SUSY Ward identities in N=l SYM theory on the lattice 
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The SUSY Ward identities (Wis) for the N=l SU(2) SUSY Yang Mills theory discretized on the lattice with 
Wilson fermions (gluinos) are considered. The study is performed in the framework of a Monte Carlo simulation 
of the model with light dynamical gluinos. The renormalization and mixing constants of the lattice SUSY current 
Zs and Zt and the additively renormalized gluino mass ms are unknown parameters of the SUSY Wis. Using 
suitable on-shell combinations of the Wis, the ratios Zt/Zs and ms/Zs are determined non-perturbatively at 
one value of the coupling constant go and two values of the hopping parameter k. 



1. INTRODUCTION 

The non-perturbative regime of SUSY gauge 
theories is highly interesting, since among other 
things it may provide a possible mechanism for 
dynamical SUSY-breaking. An increased un- 
derstanding of the non-perturbative phenomena 
of SUSY gauge theories can be achieved in the 
framework of the lattice regularization. An imme- 
diate difficulty arises however because the lattice 
regularized theory is not supersymmetric, since 
the Poincare invariance (which is a sector of the 
superalgebra) is lost. 

The present computational resources allow now 
to approach the simplest of the SUSY gauge the- 
ories, the N = 1 supersymmetric Yang Mills 
(SYM) models, where the N%-1 gluons are accom- 
panied by neutral fermionic partners (gluinos) 
in the same adjoint representation of the color 
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group. The Wilson discretization of the SYM 
models, proposed by Curci and Veneziano 0, ex- 
plicitly breaks SUSY by the Wilson term and, 
softly, by a gluino Majorana mass term. SUSY is 
only recovered upon tuning the hopping parame- 
ter to the massless gluino limit and by taking the 
continuum limit. 

An extensive analysis of the low-energy behav- 
ior of the N=l SU(2) SYM model with Wilson 
fermions was performed by this collaboration in 
a Monte Carlo simulation with light dynamical 
gluinos, see |p| and references therein. The is- 
sue of the restoration of SUSY in the massless 
gluino limit is now addressed, as proposed in j| , 
by the non-perturbative analysis of the on-shell 
lattice SUSY Wis. The same issue, from a per- 
turbative point of view, has been dealt with in 
another contribution |§. 

An approach to the SYM models with domain 
wall fermions was recently considered in [^) . 
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2. LATTICE FORMULATION 

In the Curci-Veneziano approach jl]], which we 
adopt, the standard Wilson discretization is ap- 
plied to the SYM model, where SUSY and the 
(anomalous) U(l) chiral symmetry of the contin- 
uum theory are explicitly broken by the Wilson 
term and a gluino Majorana mass term; SUSY 
is also broken by the space-time lattice. The is- 
sue of the recovering of the appropriate symmetry 
in the continuum limit is addressed through the 
discussion of the relevant Wis. 

The gluonic action is the standard plaquette 
action 



S s = /?£(l-^ReTr£/ pi ) 



(1) 



where the bare gauge coupling is given by (3 = 
2N c /gy. The fermionic action is (with Wilson 
parameter 7'=1) 



-1^ 



+ (mo + 4)A X AJ , 



(2) 



A = J2 r T r \ r being a Majorana field (gluino field) 
transforming according to the adjoint representa- 
tion of the gauge group (T r are the generators). 

2.1. SUSY WARD IDENTITIES 

The effect of the explicit breaking of the chi- 
ral symmetry by the Wilson fermions is well un- 
derstood in the framework of lattice QCD R], 
where the Wis of the continuum can be recov- 
ered through renormalization of the chiral cur- 
rents and the quark masses. The lattice SUSY 
Wis of the SYM models are derived in JIJ us- 
ing the same theoretical setup. Specifically, the 
SUSY current (x) , considered in on-shell corre- 
lation functions, undergoes multiplicative renor- 
malization and mixing with a gauge- invariant cur- 
rent T^(x); the renormalized current reads 



which, for a gauge invariant local operator 0(y) 
(y ^ x), satisfies the integrated Wis 

^(VoSo = m s J2( X (x)0(y))+0(a) 

X X 

(4) 

(the exact expressions for S fl (x), T^(x), x( x ) are 
given below). The quantity ms is obtained from 
too by additive renormalization, to 5 = too — 
ws(to ,5o)- As it is evident from the Wis (||^), 
the renormalizations Z$, Zt and rhs reabsorb the 
main lattice artifact of the explicit SUSY break- 
ing; the residual 0(a) breaking terms on the r.h.s. 
of ([!]) imply indeed effects vanishing exponen- 
tially with go ■ If these are negligible, the con- 
dition ms—0 ensures the restoration of SUSY on 
the lattice. When ms=0, also the renormalized 
gluino mass multiplying the soft breaking term in 
the lattice chiral Wis vanishes |ij as expected. 

For the numerical analysis of the SUSY Wis we 
consider two possible lattice forms of the SUSY 
and mixing currents; a point-like discretization: 

per 



T«(z) = 2VVTr[Pj$,A x ] 



(5) 



and a point-split one 
1 



2^ 

pa 



CTp CT 7 /1 Tr[p]^ ) CT [/| jll A x+( af7 x 



(*) 



1 p( ci ) 77 \ rrt 1 
' ± X +p,pa u xp n x^ X p,\ 

V 

The values of Z$ and Zt are expected to be dif- 
ferent for these two choices. Consistently with 
parity and time-reversal, P x ^ l J v is defined as a 
clover-symmctrized lattice field tensor; moreover, 
the lattice derivative in (|]) is defined, for the dis- 
cretization (|): VjL^/O) = l/2(f(x + p) - f(x - 
p,))) in the case @: V^f(x) = f(x) - f(x - ft). 

Finally, the lattice form of the operator x( x ) m 
the soft-breaking term is 



Sn(x) = Z s S,i(x) + Z T T^(x) 



(3) 



p{cl) \ 



(7) 
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Figure 1. ms/Zs as a function of the time- 
separation t at k=0.1925. 




Figure 2. Zt/Zs as a function of the time- 
separation t at k=0.1925. 



3. NUMERICAL ANALYSIS 

We use two sets of gauge configurations at 
/3=2.3 on a 12 3 x 24 lattice, at re=0.1925 and at 
k=0.194 with a statistics of 4212 and 2034 con- 
figurations, respectively. The second set is closer 
to the critical point of the chiral phase transition 
k c =0. 1955(5) as determined in All configu- 
rations were obtained by applying the two-step 
dynamical-fermion multi-bosonic algorithm, see 
|H and references therein for details. 

We consider for our study the Wis obtained by 
inserting in (Q) the two gauge invariant operators 

O^(x) = t, , Tr /', ( ,\, (8) 

i<j 

{2) {x) = ^VoiTr^Aj , (9) 

i 

with the clover and the simple-plaquette defini- 
tion of the lattice field tensor P x ^ v (the latter 
only for the operator (g)). The above operators 
undergo a (common) multiplicative renormaliza- 
tion and the form (Q) of the Wis is consequently 
not spoiled. Since spatial integration breaks hy- 
percubic symmetry in (^) these two operator in- 
sertions amount to independent choices. 

Because of the Dirac structure of the operators 
involved, Eq. (|j) gives 16 (real) Wis for every 
time-separation t = xq — yo and a given insertion 
operator O(x). Due to the Majorana nature of 



the field A T and the discrete symmetries of the 
action (0,0), only two equations are however in- 
dependent. Consequently, determinations of the 
ratios ms/Zs and Zt/Zs can be obtained, for a 
fixed 0(x) and for every t, by solving a two-by- 
two linear system. These determinations contain 
O(a) effects. 

In order to get a significant signal for the corre- 
lations, a combined APE smearing for the gluon 
field and Jacobi smearing for the gluino field was 
performed on the insertion operator. We use two 
sets of smearing parameters (set A and B) giving 
equally good results for the ratios. 

In Figs. [I] and || the determinations of the ra- 
tios ms/Zs and Zt/Zs are reported, at k=0.1925 
as a function of the time-separation t, for the two 
discretizations of the SUSY and mixing currents 
(||) @ and the two insertion operators (|8|) (||) 
(with clover field tensor). A plateau can be ob- 
served for ms/Zs (Fig. |l|) for t > 3 with the 
insertion operator 0^'(x) and for t > 4 with 
(x): in these regions of t contact terms due to 
the time-extension of the operators entering the 
correlations are absent. The signal for Zt/Zs 
(Fig. ||) is noisier. 

The results of the global fit over the range of 
time-separations t > 4 are reported in Table [jj 

Different discretizations of the SUSY or mixing 
current give different renormalizations Zs and Zt 
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Table 1 

Results for ms/Zs and Zt/Zs with the discretizations 1 and 2 of the currents and the two insertion 



operators, with plaquette and clover field tensor and smearing parameters A and B. 



K 


current 


field tensor 


smearing 


m s /Z s (0W) 


ms/Z s 


Z T /Z S (0W) 


Z T /Z S (OW) 


0.1925 


1 


plaquette 


A 


0.155(11) 




0.184(28) 




0.1925 


1 


plaquette 


B 


0.173(11) 




0.160(25) 




0.1925 


1 


clover 


A 


0.170(13) 


0.144(18) 


0.244(35) 


0.29(6) 


0.194 


1 


clover 


A 


0.118(17) 




0.190(45) 




0.1925 


2 


plaquette 


A 


0.152(9) 




-0.010(33) 




0.1925 


2 


plaquette 


B 


0.168(9) 




-0.044(32) 




0.1925 


2 


clover 


A 


0.170(12) 


0.132(16) 


-0.018(43) 


0.11(7) 


0.194 


2 


clover 


A 


0.130(15) 




-0.03(5) 





for any finite go- In particular, the discretization 
(^) gives a mixing coefficient Zt compatible with 
zero according to our data, but in case of (|[) Zt 
seems to be small and non-zero. 

Results from the two independent insertion op- 
erators OW(i) and 0^{x) differ by O(a) effects. 
The discrepancy appears to be comparable to the 
statistical uncertainty, even if the insertion opera- 
tor O^ 2 ' (x) seems to give a lower value for ms /Zs 
compared to 0^\x). Observe that also results 
coming from different lattice forms of a given op- 
erator differ by discretization errors. 

Theoretically, ms/Zs should decrease when 
the hopping parameter approaches the expected 
critical point k c , and vanish at k c ; Zt/Zs should 
take its asymptotic value for massless gluino. 
Comparison between data at k = 0.1925 and 
k = 0.194 seems to confirm this expectation. The 
statistical indetermination still prevents an accu- 
rate extrapolation of k to massless gluino. 

4. CONCLUSIONS 

The present study shows that the extraction 
of the ratios ms/Zs and Zt/Zs from the on- 
shell SUSY Ward identities is technically feasi- 
ble with the computing resources at hand. The 
main technical difficulty (related to SUSY) is that 
high-dimensional operators with a mixed gluonic- 
fermionic composition must be considered, in- 
troducing relatively large statistical fluctuations. 
This difficulty can be handled with an appropri- 
ate smearing procedure. The non-perturbative 
determination of the ratio ms/Zs can be used 



for an independent extrapolation to the criti- 
cal hopping parameter corresponding to massless 
gluinos. Our results for ms/Zs and Zt/Zs are 
consistent with the Wis (|],(|) with O(a) effects 
comparable to the statistical indetermination. 
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